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ABSTRACT 

Recently X-ray observations have shown the common presence of compact galac- 
tic coronae around intermediate-mass spheroid galaxies embedded in the intraclus- 
ter/intragroup medium (ICM). We conduct 2-D hydrodynamic simulations to study 
the quasi-steady-state properties of such coronae as the natural products of the on- 
going distributed stellar feedback semi-confined by the thermal and ram pressures of 
the ICM. We find that the temperature of a simulated corona depends primarily on 
the specific energy of the feedback, consistent with the lack of the correlation between 
the observed hot gas temperature and K-band luminosity of galaxies. The simulated 
coronae typically represent subsonic outflows, chiefly because of the semi-confinement. 
As a result, the hot gas density increases with the ICM thermal pressure. The ram 
pressure, on the other hand, chiefly affects the size and lopsidedness of the coronae. 
The density increase could lead to the compression of cool gas clouds, if present, and 
hence the formation of stars. The increase also enhances radiative cooling of the hot 
gas, which may fuel central supermassive black holes, explaining the higher frequency 
of active galactic nuclei observed in clusters than in the field. The radiation enhance- 
ment is consistent with a substantially higher surface brightness of the X-ray emission 
detected from coronae in cluster environment. The total X-ray luminosity of a corona, 
however, depends on the relative importance of the surrounding thermal and ram 
pressures. These environment dependences should at least partly explain the large 
dispersion in the observed diffuse X-ray luminosities of spheroids with similar stellar 
properties. Furthermore, we show that an outflow powered by the distributed feedback 
can naturally produce a positive radial gradient in the hot gas entropy, mimicking a 
cooling flow. 

Key words: method: numerical ISM: kinematics and dynamics X-ray: galaxies: clus- 
ters 



1 INTRODUCTION 

Though consisting of primarily old stars, galactic spheroids 
(bulges of Sb-Sa spirals as well as SO and elliptical galaxies) 
are a major source of stellar feedback in form of mass loss 
and Type la supernovae (e.g., Ciotti et al. 1991; Knapp et 
al. 1992; Mannucci et al. 2004). The specific energy of this 
feedback predicts that it should present primarily in X-ray- 
emitting hot gas. Indeed, such hot gas has been detected 
in and around spheroids, which typically contain little cool 
gas. However, it has been shown repeatedly that the X-ray- 
inferred hot gas mass and energy are far less tha n the empir- 
ical predictions from the feedback inputs (e.g.. [David et al. 
2006 ; iLi et alJliooll ; iLi and Wanell2007l : iLi et alJl2007tlWanel 
20101 ). This missing stellar feedback has most likely escaped 
into large-scale galactic halos, where the gas becomes too 
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tenuous to be detected in exi sting X-ray imaging observa- 
tions (e.g., iTang et al.l l2009al and references therein). The 
implication of this scenario is profound, because the injec- 
tion of the mass and energy into the halos could strongly 
affec t the ecosystem of the galaxies and hence their evolu- 
tion l|Tang et alj|2009gj ). 

While the above qualitative picture seems clear, there 
are key issues that still need to be addressed to understand 
both the stellar feedback itself and its interplay with the 
environment. Significant uncertainties are still present in 
the mass and energy input rates of stars (§ 2). The en- 
ergy input rate is inferred from observations of SNe in a 
large sample of galaxies of diver se optical and near-I R lu- 
minosities as well as types (e.g.. [Mannucci et al.|[2005h . as- 
suming a certain explosion mechanical energy deposited into 
the interstellar space. The stellar mass loss rate is based on 
the modeling of the 12 /im emi ssion from the circ umstel- 
lar medium of evolved stars (e.g-. lKnapp et aTlll992[ ). These 
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semi-empirical rates, uncertain by a factor of at least ~ 2 for 
individual galaxies, can in principle be directly constrained 
by the measured temperature and luminosity of galactic 
coronae. Indeed, detailed simulations have been conducted 
for relatively isolated "fi eld" spheroids and have been com- 
pared with observations (|Tang et a l. 2009a, bt lTang fc Wand 
|2010| ). which have led to a qualitative understanding of the 
feedback processes and effects on X-ray measurements. In 
particular, the feedback model expects that the specific en- 
ergy should not change substantially from one spheroid to 
another, which is consistent with little correlation between 
the measured temperatures and K-ba nd luminosities Lk 
(e.g., |Pavid et al.ll2006l ; ISun et al.ll2007l ). But the measured 
temperatures ( < 1 keV) are substantially lower than the 
expected value from simulations. The measured iron abun- 
dances of coronal gas are also typically lower than expected 
for the la SN-enriched gas (§ 2). At least part of these dis- 
crepancies can be accounted for by various 3-D effects of the 
la SN heating, which produces very low-density, hot, fast- 
moving, and enriched bubbles that hardly radiate. When the 
gas in these bubbles finally mixed with the material from the 
stellar mass loss at large radii, the X-ray emissi on becomes 
too we a k and diffuse to be e ffectively detected (|Tang et alj 
l2009al lbl; Ffang fc Wandboid ). The observed X-ray emission 
thus gives only a biased view of the coronae. We expect that 
this bias should be minimal for a corona in the ICM, how- 
ever. The high external thermal and ram pressures tightly 
confine such coronae, resulting in a low outflow speed and 
hence relatively local mixing of the feedback materials. Fur- 
thermore, one can better characterize the ICM environment 
from observations, important for a self-consistent modeling 
of a corona. Therefore, coronae embedded in the ICM are 
ideal sites to better constrain the feedback and its interplay 
with the environment. 

There have also been significant efforts in study- 
ing stellar feedback-powered coronae embedded in the 
ICM, mostly f ocusing on the ram- p ressure stripp i ng of 
hot gas (e.g., Acreman et alj 120031 ; IStevens et alj 1 19991 ; 
iToniazzo fc Sch indlcr 200j]) • It is shown that the global mor- 
phological and integrated properties, such as gas mass and 
luminosity, are strongly influen ced by the environm ent. A 
set of 2D simulations done by IStevens et al.l () 19991 ) show 
that galactic coronae can be maintained by stellar feedback 
in poor clusters whil e be efficiently stripped in rich ones. 
lAcreman et al.1 (120031 ) and lToniazzo fc Schindierl (120011 ) sim- 
ulated galaxies falling into clusters and demonstrated that a 
galactic corona reached a cyclic "stripping replenishment" 
dynamics due to the periodic orbital motion of the host 
galaxy as well as the competing processes such as stripping 
and stellar feedback. 

We focus on modeling coronae in and around 
intermediate-mass spheroids that are embedded in the ICM. 
Such a system is relatively simple with minimum effects due 
to the feedback from AGNs and to the radiative cooling of 
hot gas. We expect that the hot gas is in a quasi-steady, 
subsonic outflow semi-confined by the thermal and ram- 
pressures of the ICM. This state should be only sensitive 
to the local properties of the ICM (see § 2 for more discus- 
sion), avoiding large uncertainties in modeling the history 
of the galactic feedback and the evolution of t he environ- 
ment, as would be needed for a field spheroid l|Tang et al.l 
2009a). The simulations can also be compared with an in- 



creasing number of X-ray detection s of galactic coronae of 
such spheroids (e.g. JSun et al.l2007h . leading to an improved 
understanding of the feedback itself and its interplay with 
the environment. In particular, we examine the dependence 
of the corona properties on the specific energy of the stel- 
lar feedback and on the thermal and ram pressures of the 
ICM and check how measurements (e.g., temperature, sur- 
face brightness, overall luminosity and morphology) may be 
made to infer the parameters that cannot directly observed 
(e.g., the ICM thermal and ram pressures local to an indi- 
vidual galaxy). Here, we will present 2-D simulations only, 
which allow for an efficient exploration of a large parameter 
space. The paper is organized as follows: We briefly describe 
our numerical model and setup in § 2 and present results in 
§ 3; We discuss their implications in § 4; Finally in § 5, we 
summarize our conclusions. 



2 MODEL AND SIMULATION SETUP 
2.1 Model Galaxies 

Our model galaxy is composed of a stellar spheroid compo- 
nent and a dar k matter halo. W e use the spherical Hernquist 
density profile (|Hernquistl 19901 ) to represent the stellar mass 
distribution: 
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where M s is the total stellar mass, and a is the scale radius. 
This density profile results in a gravitational potential 

^) = -tS £ T- ( 2 ) 
[r + a) 

The above stellar mass distribution resembles the de Vau- 
couleur's Law; The relation between the half-light radius R e 
and the scale radius is R e — 1.8513a. 

We characterize the d ark matter halo with the NFW 
profile (|Navarro et al.lll997l ). 
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where r d is the scale radius of the dark halo, and po is defined 

as 
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in which p cr it is the critical density of the universe. The dark 
halo has a mass M vir within the virial radius r vir , which is 
defined to have a density that is A V i r times the mean matter 
density of the universe p C r«f2 m . We adopt A V i r to be 250. 
Therefore, we have the relation 

47T 3 

Ml v i r — — T v i r A v i r Pcrit^m , (5) 

where c is the concentration factor defined as c = ^-f^- and 
is related to M vir . We set c = 13 ijEke et al.ll200ll ). Thus, 
for a given cosmology and a given M v i r , the mass profile of 
the dark halo is totally determined. 



2.2 Stellar Feedback 

The stellar mass and energy feedback in spheroids are dom- 
inated by the mass loss and la SNe of evolved stars, respec- 
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tively. We neglect the energy input from the random motion 
of stars and hence their ejecta. The total energy released by 
type la SNe is 



E = E SN n SN (^f- 



(6) 



where ngjy = . 00035 yr~ for E/SO galaxies according to 
IMannucci et al.l l|2005l ). The empirical mass input rate from 
the stellar mass loss is 



M = ml 



where m = 0.021 M Q yr 1 according to lKnapp et al. 
Assuming that the mechanical energy of each SN is 10 51 erg, 
the specific energy of stellar feedback is /3 = S- ~ 5 keV per 
particle. To account for the uncertainties in these rates and 
assumptions, we also sample three different lower values of 
the specfic energy for comparison with observations (Table 
1). We assume that the energy and mass inputs follow the 
distribution of the stellar mass. 

In addition, each la SN produces ~ 0.7Mq of iron 
ejecta. We assume that the iron abundance of the mass loss 
from stars is solar. If the ejecta is fully and instantaneously 
mixed with the mass loss, the expected iron abundance rel- 
ative to the solar value is then 
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However, in observa tions, supersolar metallicity is quite rare. 
iTang et alJ (|2009tJ ) and ITang fc Wanel 1)20101) have shown 
that la SN ejecta may not be efficiently mixed with stellar 
mass loss on microscopic scale, resulting in a low effective 
metallicity of the ISM. While our focus is on the environ- 
mental effect. In our simulation, we set the iron abundance 
of the input mass to the solar value. 



2.3 Simulation Setup 

The gas hydrodynamics and metal abundance distribution 
with both the stellar feedback and radiative cooling can be 
described by the following equations: 



dt 



+ V • 9 v) 



d(p 9 v) 
dt 

d(p g e) 
dt 

dt 



+ V-(p 9 vv + P[/]) = 
V-[(p 9 e + P)v] = 

+ V ■ (pgXironV) = 



S n 



S e - n e n l A(T, Z) + p 9 g ■ v 

Siron- (9) 



The first equation is mass conservation law, with p g denot- 
ing the mass density of the coronal gas. The second is mo- 
mentum equation. P is gas pressure and g is gravitational 
acceleration. The third one is energy equation, e stands for 
the specific energy of the gas, including both thermal and 
kinetic components. The second term to the right is cooling 
rate. W e adopt the cooling curve from lSutherland fc Dopital 
(1993), assuming an optically-thin thermal plasma in colli- 
sional ionization equilibrium. For calculation of the radia- 
tion in a specific band, we use the Mekal model, extracted 
from the X-ray spectral analysis software XSPEC. We use 
the fourth equation to keep track on the iron mass fraction, 
which is denoted by Xiron- 



We conduct ou r simulations with the FLASH code 
l|Frvxell et al l [2000), an Eulerian astrophysical hydrody- 
namics code with the adaptive mesh refinement (AMR) ca- 
pability. The simulated region is fixed in the galaxy-rest 
frame using cylindrical coordinates, with z ranging from 
— 50 kpc to 50 kpc and the radius from to 50 kpc. The axis 
of the cylinder is through the center of a simulated spheroid 
and along the direction of its motion. The upper and lower 
boundary conditions are fixed so that the ICM flows in and 
out the simulation region at a constant speed. This mimics 
the motion of galaxies through a local cluster environment. 
We apply reflection boundary condition at r = and diode 
boundary condition, which only allows gas to flow out, to 
the right side of the simulation region. Compared with the 
simulation box, which is 50 kpc by 100 kpc, the coronae are 
only on the order of 1 kpc to 10 kpc across. Data near the 
outer regions will be excluded in our analysis to avoid any 
potential artifacts introduced by the assumed outer bound- 
ary condition of the simulations. Also, the outer region will 
not be shown in the following images. We allow the resolu- 
tion to reach 0.1 kpc, so that the small coronae can be well 
resolved. 

When a simulation starts, there is no interstellar gas in 
the galaxy. As the simulation progresses, the stellar feedback 
gradually accumulates in and around the spheroid to form a 
corona, which is characterized by its higher iron abundance. 
In the mean time, the ram-pressure and turbulent motion 
strips gas at the outer boundaries of the corona. We end the 
simulation when it reaches a statistically quasi-steady state. 
Empirically, the time to reach such a state is r 9 ~ 0.2 Gyr, 
while the time for the ICM to pass the simulation region 
ranges from 0.05 Gyr to 0.2 Gyr. Representative results are 
all extracted from the simulations after this time. 

The presence of a local quasi-steady state is a reason- 
able assumption for a compact corona. As a galaxy moves 
through a cluster, the ICM condition can of course change 
drastically. But the time scale for such a change is typically 
much longer than the dynamic time for the corona to adapt 
the local environment. For a cluster of a characteristic size of 
~ 1 Mpc and temperature of 2 keV, for example, the cross- 
ing time for a galaxy moving roughly at the sound speed is 
t c i ~ 2 Gyr. In contrast, for a corona of a typical size ~ 5 
kpc and temperature ~ 0.8 keV, the sound crossing time 
is only r C 2 ~ 20 Myr. Even if a corona is totally destroyed 
at some point(e.g. at the central region of a cluster), the 
re- building time scale r 9 , as dicussed above, is still shorter 
than the environment change time scale. Therefore, the local 
quasi-steady state is a reasonable assumption for character- 
izating the environmental impact on galactic coronae. 



3 RESULTS 

We have simulated a set of cases to characterize the depen- 
dence on key parameters. Table 1 lists our adopted model 
parameter values. The different combinations of the ICM 
density (rii, the number density of all particles) and temper- 
ature (Tj) as well as the Mach number (Mk) an( f specific 
energy (j3i) of the model galaxy form a set of 48 cases. 

Here we present the gas properties extracted from the 
simulations. We first detail the results for a representative 
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Table 1. Model Parameters 



Model galaxy 


Stellar Mass(10 11 M Q ) 


2.0 




Dark Halo Mass (10 11 M ) 


10 


ICM properties 


ICM density (1CT 4 cm" 3 ) 


3.3(m), 10(n 2 ) 




ICM temperature (10 7 K) 


2.0(Ti), 6.0(T 2 ) 




Iron abundance ( Zq) 


0.3 




Mach number 


0.6(Mi), \.2{Mi), 1.8(A4 3 ) 


Stellar feedback 


Mass loss rate ( M Q /10 n M Q /yr) 


0.32 




Specific energy (keV) 


1.2(/3i), 1.8(/3 2 ), 3.0(/3 3 ), 4.8(/3 4 ) 




Iron abundance ( Zq) 


1.0 



case n\TiM2f$2 (§ I3.1|l and then discuss the similarity and 
significant variation among the different cases (§ 13. 2p . 

3.1 Representative Case 

Fig.[T]shows a snapshot of the representative simulation case 
niT2A / f2/?2 in terms of the Mach number, thermal pressure, 
density, and temperature distributions. At the very front of 
the corona is a smooth and distinct boundary that sepa- 
rates the corona from the ambient medium. This is a con- 
tact discontinuity, across which the density, temperature and 
metallicity change abruptly. Compared to the surrounding 
medium, the corona is cooler and denser. The iron abun- 
dance inside the corona is a constant which is equal to the 
value of injected material (Fig. [2}. Outside the corona the 
abundance drops rapidly to the value of the ICM, although 
it is contaminated by the local stellar feedback. Therefore, 
we can use the abundance to trace the morphology of the 
corona gas. Inside the corona, the Mach number of the gas 
flow is so low (~ 0.1) that it is almost hydrostatic (see § 14.11) . 

While the main body of a corona can reach a nearly 
steady state, both the interface with the surrounding ICM 
and the tail are unstable. The individual features in these 
later parts can strongly fluctuate with time. The side horns 
are characteristic sign of the Kelvin-Helmholtz (KH) insta- 
bility. As a result, the corona gas is torn off and pushed 
back to form a chaotic tail. Therefore, the stripping is pri- 
marily due to the hydrodynamic instability rather than the 
ram-pressure itself. Similarly, the instability also leads to 
the dynamic mixing of the corona gas with the ICM, al- 
though numerically this is achieved on the spacial scale of 
the simulation resolution. 

Fig.[3]shows the 1-D distribution of the iron abundance, 
density, temperature and entropy of the corona along the z 
axis of the simulation box. Here the entropy is defined as 
5* = n -^-i (7 is the specific heat ratio of the gas and is equal 
to |). The distributions are averaged over a time span of 50 
Myr when the simulation has reached a quasi-steady state. 

While the gas density drops substantially from the cen- 
ter of the spheroid to the outer outskirt of the corona, the 
temperature does not change much (Fig. [3]). The specific en- 
ergy of the feedback determines the specific enthalpy of the 
corona gas and therefore the temperature (T = 2 A ) when 
the radiative cooling is not important as in the present case. 
The small drop of the temperature towards the outskirt (by 
a factor of up to 1.2) is largely due to the outflow that needs 
to climb out of the gravity potential. But, because of the dis- 
tributed nature of the mass and energy injection, the drop 
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Figure 1. Model n\TiMifi'z seen in the Mach number, thermal 
pressure, density, and temperature. 



is much smaller than what is predicted from the Bernoulli's 
law for ideal gas moving from the center to the outskirt. 

Another interesting characteristic of the simulated 
corona is the positive radial entropy gradient. This is ap- 
parently caused by the nearly constant temperature profile 
and the steep density drop from the center to the outskirts, 
as required to maintain a nearly hydrostatic state of the 
corona (§ 14.10 . Physically this positive entropy gradient is a 
natural result of an outflow that is continuously heated by 
the stars along the way out. 

3.2 Similarity and Variance among the Cases 

Here we focus on the similarity and significant variance in 
the hot gas properties of the various simulated cases, in ref- 
erence to the representative one (niT 2 A / (2/3 2 ) detailed above. 
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Figure 2. Iron abundance map for model nxTiMifli- 



In most of the simulation cases, the coronae are clearly 
in the outflow state. The radiative cooling is not important, 
except for some such as niT\M\fi\, ri2T\M.%fi\, niTzMiPi, 
nzTzMxfix and niTiMifix , with combination of low spe- 
cific energy and high ICM pressure. The gas at the spheroid 
center is so dense that an inward cooling flow is developed 
in the inner region. (Fig. 2}. Such cooling flows, commonly 
seen in similar models and simulations, may naturally in- 
duce activities of the central supermassive black holes. The 
feedback from such activities has b een proposed to sub- 
stantially reduce the net cool ing (e.g. iMathews fe BrighentH 
l2005l ; lFabian fc Sandersj|2009l ). While a study of this topic is 
beyond the scope of the present work, we here keep our fo- 
cus on discussing the morphological and physical properties 
of the outflow cases, which probably represent more typical 
cases of galactic coronae in the ICM. 

For all the 712T2.M3 cases, coronae fail to form due to 
the high thermal and ram pressures, leaving only a track 
with high iron abundance (Fig. [5}. In each of these cases, 
the density peaks away from the spheroid center. 

Fig. [6] show outlines of the coronae. These outlines 
are represented by the iso-abundance contours of the value 
of Zq, or approximately the contact discontinuity between 
corona gas purely ejected by the spheriods and the ICM. 
The characteristic size of a corona is sensitive to the ther- 
mal pressure of the surrounding ICM as can be seen in panel 
(a) of Fig. [6] The corona in the ICM of n\T\, which is typical 
of a group or poor cluster, is more than lOkpc across, while 
the corona in the ICM of 712T2, which is typical of the core 
region of a relatively rich cluster, is less than 3kpc across. 
The motion of the host galaxy relative to the surrounding 
medium, which is represented by the Mach number, mainly 
influences lopsidedness of a corona. In the subsonic case, the 
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Figure 3. Iron abundance, denisty, temperature, and entropy 
distributions along z-axis of the simulation of case n\T2-Mif^2- 
The thick solid lines are time-averaged profiles over a time span 
from 250 to 300 Myr (the residual wigglers are present due the 
the limited number of snapshots used in averaging), while the 
thin lines are extracted from several snapshots. The horizontal 
red line in the temperature panel corresponds to T = 2 ^ k . 
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Figure 4. Iron abundance and density distributions in the 
niT2MiPi case. The arrows represent the velocity field. Note 
the central density peak (the red spot) at the very center. 

corona is almost spherical, while in the supersonic case with 
the Mach number as high as 1.8, the corona is significantly 
narrowed and elongated. As the input energy of the stellar 
feedback increases, the lopsidedness increases. This is not 
surprising. As the stellar feedback becomes more energetic, 
the density of the corona gas will be lower. As a result, the 
gravitational restoring force become less important, com- 
pared with the ram pressure. 
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Figure 5. Iron abundance and density distributions in the 
rtiT2.Mi/3i case, in which the corona fails to form. 
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Figure 7. Temperature profiles along z-axis. The red horizontal 
lines in (c) represent 2~5k^' '^ le re f erence model nxTi-Mifil is 
represented with solid line in all the three figures. Panels (a), 
(b) and (c) correspond to variation in nT, Mach number and f} 
respectively, illustrating how the coronae respond to the changes 
of the surrounding environment and stellar feedback. 



10 



5 



-5 



(a)' 


_ 'n,T,^ 2 M 2 


(bj 


_ _ n^p^M, 


(c) 


' 'n 1 T 2 ^ 2 M 3 




_ _ n 2 T,p 2 M 2 




n,T 2 p 2 M 2 




- - n iT 2 p 3 M 2 




n,T 2 p 2 M 2 




n,T 2 p 2 M 3 




__.nJ 2 p 4 M 2 




_ _ n 2 T 2 p 2 M 2 


























- j 1 







02468 10 2468 10 2468 10 
r /kpc p /kpc p /kpc 

Figure 6. Outlines of the simulated coronae. The outline for each 
case is defined to be the contour of the iron abundance equal to 
Zq. The reference model n\T2hA2&2 is represented with solid 
line in all the three figures. Panels (a), (b) and (c) correspond 
to variation in nT, Mach number and /3, respectively, illustrat- 
ing how the coronae respond to the changes of the surrounding 
environment and stellar feedback. 

The dependence of the corona temperature on the spe- 
cific energy of injected material can be clearly seen in Fig. [7] 
The peak value of the corona temperature is always roughly 
2 j? k ; it drops slightly outwards largely due to the presence 
of gravitational potential as described in § I3.ll 

We plot the peak density values of the coronae in Fig. [8] 
and selected density profile along the z-axis in Fig. [9] illus- 
trating how the coronae respond to the changes of surround- 
ing environment and inner stellar feedback. Generally, the 
corona density depends strongly on the thermal pressure of 
the surrounding medium, but only weakly on the ram pres- 
sure. In our simulation, the the Mach number ranges from 
0.6 to 1.8, with the ratio between the largest ram pressure 
and the smallest one as high as 9, the corona density changes 
slightly except for case T^TbA^iA an d its high Mach num- 
ber version n2T2M.-2.fi1- This makes sense because unlike 
thermal pressure, which compresses a corona from all direc- 
tions, ram-pressure only acts on the front side and therefore 
mostly pushes the gas backwards rather than compressing 
it. The coronae in cases of 712 T\ are compressed more than 
those in cases of niTi, although the ICM thermal pressure 
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Figure 8. Peak density values of the coronae. Cases with central 
cooling flow and those in which the coronae fail to form are not 
shown. Cases with different ICM thermal states are represented 
with different symbols, while different Mach numbers are coded 
in different colors with red for 1.8, green for 1.2 and blue for 0.6. 
To avoid overlap among the symbols, the higher Mach number 
models are shifted to the right a litte bit and the subsonic cases 
to the left. 



is the same. This is caused by larger focusing effect of grav- 
itational force of the galaxy on the ICM with lower temper- 
ature. 

Fig.[l0]includes the luminosity of each simulated corona 
in the 0.3 — 2.0 keV X-ray band. For simplicity, we use the 
iron abundance distribution to define the shape of a corona 
as what we do to plot the outline of a corona. Every grid 
point with iron abundance equal to the solar value is in- 
cluded to measure the luminosity. The emissivity is a func- 
tion that depends on both temperature and metallicity. If 
we assume the metals, mostly iron, produced by la SNe, are 
fully mixed in the corona gas, the total luminosity would be 
enhanced by a factor of about 3. It is clear that the lumi- 
nosity decreases with the increasing Mach number and/or /3. 
But the dependence on the ICM thermal pressure is not that 
simple: the luminosity tends to increase with the pressure in 
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Table 2. a for various cases. 
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Figure 9. Density profiles along z axis. The reference model 
rtiT2.M2/32 is represented with the solid lines. Panels (a), (b) 
and (c) show the dependence on nT, Mach number, and /3, re- 
spectively. 
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Figure 10. Luminosity of the coronae. Cases with central cooling 
flow and those in which the coronae fail to form are not shown. 
Cases with different ICM thermal states are represented with dif- 
ferent symbols, while different Mach numbers are coded in differ- 
ent colors with red for 1.8, green for 1.2 and blue for 0.6. To avoid 
overlap among the symbols, the higher Mach number models are 
shifted to the right a litte bit and the subsonic cases to the left. 



subsonic cases (lower Mach numbers), while the trend goes 
in the opposite direction in the supersonic cases. 



4 DISCUSSION 

The above results give a basic characterization of galactic 
coronae powered by stellar feedback and semi-confined by 
the thermal and ram prssures of the ICM. In this section, 
we first give a physical account of the apparent ICM im- 
pacts on the corona characteristics as described above, then 
compare the results with observations to constrain the stel- 
lar feedback, and finally discuss the implications for other 
galactic properties. 
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4.1 ICM Impacts on the Coronae 

The properties of a corona are affected by several competing 
processes, the stellar feedback, the galaxy gravitational at- 
traction, and the thermal/ram pressures of the ICM. In the 
simplest case, when the gravity and external pressure can 
be neglected, the density at the center would then be deter- 
mined entirely by the mass/energy injection. In this case, 
we can define a characteristic density as 



3 M 



Pel 



(10) 



47T a 3 

Here r = — is the dynamic time scale of the corona, where 
c s is the sound speed. Considering that the temperature is 
determined by the specific energy /3, the characterisitic value 
of the pressure is 



Pel = 



3 M a/(7-1)£ 



47T a 2 



(11) 



In the other extreme, when the gravity is important 
(i.e., the corona is nearly hydrostatic), the Mach number of 
the outflow must be low. Clearly, in this case, the external 
thermal pressure becomes important as well. Because the 
corona is nearly isothermal, the pressure distribution is 

fim p ($ - $o) 



m(P) = ln(Po) 



k B T 



(12) 



where Po is the pressure at the outer boundary and $ — <E>o 
is the gravitational potential difference. Assuming that the 
corona size is considerably larger than the scale a of the 
stellar spheroid, we obtain a characteristic peak pressure as 



Pc,2 = Picm exp 



( GMs 

V a 



7 



(7-1)0 



(13) 



To characterize the relative importance of the feedback to 
the gravity, we define a dimensionless parameter as 

_ PoA 

Pc,2 ' 



(14) 



The values of this parameter for the simulated cases are 
listed in Table. 2. Fig. [11] compares the hydrostatic solu- 
tions (Eq. I13p to the simulations with three different a pa- 
rameters. In those cases with small a, such as n2T2A4i/32, 
in which the gravity dominates over the feedback, the hy- 
drostatic solutions give nearly perfect matchs to the simu- 
lated pressure profiles. While in a case like n\T\M\fi/± (with 
a = 3.3), where the feedback dominates, the deviation of the 
hydrostatic solution from the simulated profile is apparent. 
For the cases which we think are plausible (/3=1.8keV, see 
§ I4.2[I . a is about 1.0 or much less than 1.0, and the coronae 
are largely in hydrostatic state. 

Although the same model galaxy (in terms of the stellar 
and dark matter masses) is adopted for all the simulation 
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Figure 11. Comparison between numerical pressure profiles and 
the theoretical profiles with isothermal assumption. The solid 
lines are the pressure profiles extracted from the simulations and 
the dashed lines represent the corresponding hydrostatic solu- 
tions. The region inside the coronae is colored in red. In each of 
the panels, the peaks of both of the numerical profile and theo- 
retical one are positioned at the same point. 



cases, the resultant X-ray luminosities can still differ by up 
to two orders of magnitude, due to the different choices of 
the Mach number, j3, and/or thermal pressure values of the 
ICM. Because the corona temperature is determined by /3, 
the gas density is 



Picm 



(15) 



if the corona is in a nearly hydrostatic state. Therefore, the 
surface brightness of a corona provides a measure of the 
ambient ICM pressure and thus may be used to estimate 
the line-of-sight position in a cluster. 

However, the X-ray luminosity of a corona depends 
on several factors. Fig. [TO] shows a clear anti-correlation 
between the luminosity and specific energy j3. This anti- 
correlation is primarily due to the density decrease with the 
increase of /?, although it does not strongly affect the size 
of a corona. The ICM thermal pressure tends to squeeze 
the corona, hence enhance its luminosity (subsonic cases). 
But this effect is complicated by the presence of the ram- 
pressure. As the Mach number increases, the lumonosity 
can decrease because the ram-pressure stripping reduces the 
overall size ofthe corona. These dependences on the environ- 
ment as well as the stellar feedback energetics may naturally 
explain the observed large dispersion of L x /Lk for spheroids 
of similar Lk- The complications in the dependences may 
also account for the lack of a clear observed trend in the 
ICM environment effect on X-ray luminos ities of coronae 
i|Sun et al.ll2007l ; iMulchaev fc Jeltemdl2O10h . 



4.2 Implication for the Feedback Model 

In our model of the galactic coronae, the gas temperature 
is primarily determined by the specific energy of the stellar 
feedback and thus should not change significantly with the 
stellar mass. This independence on the mass or Lk is consis- 
tent with the temperature m easurements of the coronae of 
intermediate-mass sphero i ds (IDavid et aL 2006; Sun et al.l 
120071 ; iJeltema et all 120081 ; iBoroson et al.ll2010T ). This is in 
contrast to the correlation between the temperature and 
Lb for mo re massive systems such a s clusters and groups 
of galaxies l|Helsdon fc Ponmanll2003r ). In particular, galaxy 
clusters show a well-defined scaling law between the temper- 
ature and luminosity of the observed hot gas, which is a nat- 



ural result of the predominant gravitational heating in the 
self-similar cluster formation. The scaling law for lower mass 
systems (e.g., groups of galaxies) is known to be slightly 
different from that for clusters (e.g., showing an 'entropy 
floor'), which is believed to be an imprint of preheating(e.g., 
starburst and early AGNs). Correlation between Lb and the 
temperature of hot gas is observed in massive X-ray-bright 
elliptical galaxies, especially for central galaxies in groups 
and clusters fe.g.. lO'Sullivan et al.ll2003l ). But the entropies 
are found lie below the entropy floor(~ 10 9 K cm -2 ) discov- 
ered in groups of galaxies. Radiative cooling could account 
for the low entropy, although how this runaway process may 
be balanced by the heating due to the mechanical inputs 
from both stellar and AGN feedbacks remains unclear. We 
have shown that the low value and the positive radial gra- 
dient of the entropy are expected from the distributed feed- 
back in intermediate-mass spheroids, in which the radiative 
cooling is not important. Therefore, we may conclude that 
the coronae of intermediate-mass spheroids represent the ex- 
treme case in which the stellar feedback plays a dominant 
role, which means they are produced by stellar mass loss 
and heated by SNe. 

We may further constrain the stellar feedback based 
on the measured temperature of the coronae. Though with 
a relative large dispersion, the measured temperatures are 
mostly fall in the range of 0.5 — 1.0 keV, which is signifi- 
cantly higher than those measured in field spheroids, but is 
still substantially lower than what is inferred from our 2- 
D simulations if the canonical specific energy value of the 
stellar feedback is assumed (~ 5 keV; see §0. Part of this 
discrepancy could still be due to the 3-D effects of discrete 
he ating by la SNe, a s mentioned in 8 1 (a nd characterized 
in iTang et aTll2009bl ; pTang fc Wangj feoich . But we expect 
that such effect should be substantially weaker in the com- 
pact coronae embedded in the high pressure ICM and that 
the measured temperature should more faithfully reflect the 
specific energy of the feedback. To match the measured tem- 
perature range of the coronae requires a specfic energy of 
~ 1.5 — 3 keV, or a factor of ~ 2 — 3 lower than the canon- 
ical value (5 keV; §[2|. This factor is probably still within 
the uncertainties of the semi-emipirical mass and energy in- 
jection rates. Further, the assumed mechanical energy per 
la SN could be somewhat (e.g., a factor of ~ 2) less than 
10 51 ergs. Also a considerable fraction of the energy can be 
used to generate cosmic rays, magnetic field and turbulent 
motion. The diversion of the energy into these various forms 
could significantly reduces the temperature, although the 
hydrodynamics of the coronae, hence the density and pres- 
sure distributions, should not be significantly affected. The 
simulated coronae with /3 = 1.8 keV generally have individ- 
ual luminosties of a few times 10 39 erg s _1 , consistent with 



the observed range of 10 



10 40 erg s _1 dSun et alj|2007h . 
A considerably large value of f) is not favored, because it 
decreases the expected luminosity steeply (Fig. 1 10 p . 



4.3 Implifications for Understanding Other 
Galactic Components 

We discuss here the potential impacts of the pressure or 
density enhancement of the coronae on the fueling of the 
central SMBHs and the evolution of cool gas, if present in 
the spheroids. 
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Figure 12. Estimated AGN power. Cases with cooling flow and 
those in which the coronae fail to form are not shown. Different 
ICM thermal states are represented with different symbols and 
different Mach numbers are coded in different colors with red for 
1.8, green for 1.2 and blue for 0.6. To avoid overlap among the 
symbols, the higher Mach number models are shifted to the right 
a litte bit and the subsonic cases to the left. 



The simulation shows that the central density of a 
corona is sensitive to the thermal pressure of the surround- 
ing medium. To infer the power o f a SMBH, we adopt the 
Bondi accretion rate ( Edgar 2004): 



M = 



A-kG 2 MI hP 



(16) 



where Mbh is the mass of the SMBH, while p and c s are the 
density and sound speed at the center of a corona. Assuming 
the fraction of the accretion energy released is rj = 0.1, the 
power of the SMBH can be approximated as 



,7.0 x 10 7 K 



(17) 



The SMBH mass can be estimated from its correlation with 
the s pheroid mass Mbh ~ 0.006Mb u i ge l)Magorrian et al.l 
1998). Fig. 1121 shows the dependence of the power on the 
ICM state, the Mach number of spheroid, and the specific 
energy of the feedback. 

In particular the luminosity for the most plausible 02 
cases is about 10 43 erg s _1 , which falls in the range of the 
X-ray power of a low luminosity AGN. This indicates that 
a compact corona built up by stellar feedback and embed- 
ded in cluster environment could feed a moderate AGN. The 
ICM pressure also tends to enhance the accretion, consistent 
with the finding that the galaxies with L x > 10 42 erg s _1 
AGNs are more cen trally concentrated than ones without 
|Martini et al.ll2007h . These galaxies with AGNs are not 
dominated by galaxies that have recently entered the clus- 
ters. Similar conclusions are also rea c hed i n more recent 
studies, such as the one bv lHart et all (|2009l ). based on the 
analysis of a sample of P\aghz > 3 x 10 23 WHz -1 radio 
galaxies and I/o.3-8feeV > 10 42 erg s _1 point sources. 

The pressure enhancement could also have signifi- 
cant impacts on cool gas in a galaxy. Under high pres- 
sure, cool gas exists preferentially in molecular form rather 



than atomic one. The compression of cold gas because 
of th e ICM pressure coul d further lead to star forma- 
tion (Bck ki fc Couch! 120031 ). depriving the galaxy of the 
gas further. Thus it is expected that galaxies contain less 
amounts of cool gas in clusters than in the field. These 
impacts should affect cool gas not only in spheroids, but 
in spirals as well, con sistent with existing observations 
jYoung fc Scovillelll99lh . 



5 SUMMARY 

We have conducted a range of 2-D hydrodynamic simual- 
tions of galactic coronae that result from gradual energy 
and mass feedback in stellar spheroids moving in the ICM 
enviornment. We have focused on spheroids that are in the 
intermediate-mass range (corresponding to Lk ~ 10 11 — 
10 12 Lk,o) so that both the AGN feedback and the radia- 
tive cooling of the hot gas could largely be neglected. We 
explore the dependence of corona properties on the specific 
energy of the stellar feedback as well as on the ram and ther- 
mal pressures of the ICM. Our major results and conclusions 
are as follows: 

(i) X-ray coronae embedded in clusters could be naturally 
explained by the subsonic outflow driven by the stellar feed- 
back, semi-confined by the ram-pressure and compressed by 
thermal pressure of the surrounding ICM. The corona tem- 
perature depends primarily on the specific energy of the in- 
put material in such a way that T ~ 2 ^ k . The decrease 
of the thermal energy due to the climbing of the gravita- 
tional potential and the expansion is largely compensated 
by the distributed heating by la SNe. This result naturally 
explains the lack of the correlation between the temperature 
and K-band luminosity for the spheroids in our considered 
mass range. An outflow powered by a distributed feedback 
also has a positive radial entropy profile, mimicking what 
may be produced by a "cooling flow". 

(ii) The coronal gas is typically in an approximate hydro- 
static state. As a result, the density of the corona gas de- 
pends strongly on the thermal pressure of the ICM, but only 
weakly on the ram pressure. Therefore, the surface bright- 
ness of X-ray emission is a good measurement of the thermal 
ICM pressure, which may be used to estimate the line-of- 
sight location of a spheroid in a cluster. The total X-ray 
luminosity of a corona decreases with the increase of the 
feedback specific energy. The thermal pressure tends to in- 
crease (or reduce) the luminosity in subsonic (supersonic) 
cases. 

(iii) The semi-confinement of the coronae by the ICM al- 
lows a good constraint on the energetics of the stellar feed- 
back. To be consistent with the observed X-ray luminos- 
ity and temperature, the specific energy of the feedback 
should be ~ 1.5 — 3 keV, a factor of 2-3 smaller than the 
value inferred from the commonly accepted semi-empirical 
la SN and mass-loss rates, assuming the mechanical energy 
of 10 51 ergs s _1 per SN. 

(iv) The relatively high pressure of the coronae in the 
ICM may have important implications for understanding the 
AGN activity as well as the cool gas properties in spheroids. 
The density increase caused by the ICM pressure, for exam- 
ple, could enhance the Bondi accretion, which may explain 
the observed central concentration of AGNs in clusters. The 
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high pressure can further compress the gas to form molec- 
ular clouds and enhance star formation. The combination 
of the enhanced consumption and the ram-pressure striping 
can naturally lead to the deprivation of gas and subsequent 
passive evolution of galaxies in clusters. 
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